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^ '. Abstract 



We present a systematic study of transverse mass (my) spectra of mesons and baryons at RHIC 
and LHC energies. In an earlier study, it was shown that all the mesons produced in p+p and d+Au 
collisions at ^/snn = 200 GeV follow mj- scaling while in Au+Au collisions at same energy the 



, mesons with strange and charm quark contents do not follow my scaling which can be attributed 

to medium modifications. We extend this study for baryons produced in all above colliding systems 
at ^/snn = 200 GeV. Although all available baryon spectra behave differently from mesons but 



are found to scale with protons for both p+p and d+Au collisions. In case of Au+Au collisions, 
the strange baryons behave differently from protons. This study has also been performed on the 



> 
oo 

' meson spectra in p+p collisions at 62.4 GeV which lies in between the highest RHIC energy and 



SPS energy where the my scaling was first observed and we arrive at same conclusions as those 
at 200 GeV. We test the ra-r scaling at LHC with the first meson and baryon spectra measured in 
p+p collisions at 900 GeV and find that the m? spectra of kaons and <ft do not scale with pions 



that well. Here we use all data measured in the mid rapidity region. Proton to pion ratio is studied 
as a function of mj, at different energies and for different systems and is found that this ratio for 
p+p at 900 GeV is quite low from that in p+p collisions at RHIC energies. 
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I. INTRODUCTION 



Heavy ion collisions at relativistic energies are performed to study the properties of matter 
at high temperature where a phase transition to Quark Gluon Plasma (QGP) is expected. In 
Au+Au collisions at ^/snn = 200 GeV at Relativistic Heavy Ion Collider (RHIC), many sig- 
nals point to the formation of Quark Gluon Plasma (QGP) jl-4|. RHIC continues to study 
the detailed properties of the strongly interacting matter using p+p, d+Au and Au+Au 
systems at various colliding energies from 200 GeV down up to 7.7 GeV. First results from 
Pb+Pb collisions at LHC already giving possible glimpses of QGP from precise measure- 
ments of observables such as jets p and quarkonia jf^O]- First measured charged particle 
spectra from LHC [8j are also available which give the global features of the collisions. 

Measurements of transverse momentum spectra for particles emerging from p+p collisions 
are used as a baseline to which similar measurements from heavy ion collisions are compared 
to study collective effects. Long ago WA80 collaboration found that the spectral shapes of 
Ti and r] mesons measured in S + S (200 AGeV) collisions are identical when plotted as a 
function of ttit [9]. This property is known as tjit scaling and has been extremely useful to 
extrapolate the meson spectra to unkown regions and to obtain the mesons spectra which are 
not measured. One application of such analysis is to obtain cocktail of decay products from 
produced hadrons {4, 10]. Obtaining the integrated yields of hadrons is another important 



application where uit scaling can be used. 



In an earlier study [ill ], we found that all mesons follow rrijn scaling in p+p and d+Au 
collision at ^/snn = 200 GeV, while the mesons with strange and charm quark contents do 
not follow it in Au+Au collisions at the same energy which can be associated to medium 
modifications. The more detailed centrality analysis of Au+Au collision reveal that in case 
of 77, the data of all centralities are very well reproduced by scaled pion data. The ratios 
obtained also remain approximately same for all centralities. For particles, kaon and (ft, 
peripheral collision data is better reproduced as compared to central collision data and the 
ratios also increase as the collision becomes more central. We extend this study for baryons 
produced in all above colliding systems at ^snn = 200 GeV. We study proton and strange 
baryons such as A, 5, Q and their anti-particles. Although all available baryon spectra 
behave differently from mesons but are found to scale with protons for both p+p and d+Au 
collisions. In case of Au+Au collisions, the strange baryons behave differently from protons. 



This study has also been performed on the meson spectra in p+p collisions at 62.4 GeV. 
This lies in between the SPS energy and highest RHIC energy. The conclusions at 200 GeV 
hold at 62.4 GeV. Finally and most importantly we test the ra^ scaling at LHC with the 
first meson and baryon spectra measured in p+p collisions at 900 GeV. It is found that the 
kaons and the <j) do not scale that well with pions; but in case of baryon, the ra^ scaled 
curve is in good agreement with the A and with the small data set of S. We also study the 
proton to pion ratio in different systems as a function of energy. It is found that this ratio 
for p+p at 900 GeV is quite low from that in p+p collisions at RHIC energies. 



II. FIT PROCEDURE USING m T SCALING 

In this section, we describe the procedure of fitting of hadron spectra using tjit scaling. 
To give a good descritption of pion spectra in wide range, the PHENIX collaboration 
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14l | used a simple form referred as the modified Hagedorn 



ormula. This formula has 



been extensively used for p+p, d+Au and Au+Au collisions ll| which in terms of % is 
given by: 

d 3 N A 
E- 



dp 3 exp(-am T - bm 2 T ) + ^ 

= U/p (y/Pr + m l/ P ) , (!) 

which is close to an exponential form at low pt and a pure power law form at high p?. Here 
f n / p is the pion or proton fit function, m 7T / p is the rest mass of pion or proton and A, a, b, 
p and n are the fit parameters. Then we obtain the spectra of mesons (or baryons) using 
pion (or proton) fit function as: 

J 3 N 
dp 3 



E^r = S U /p Up 2 T + ml), (2) 



where rrih is the rest mass of the corresponding meson (or baryon) . The factor S is the relative 
normalization of the meson (or baryon) rtvr spectrum to the pion (or proton) tut spectrum 
which we obtain by fitting the experimentally measured meson (or baryon) spectrum. We 
use pion fit function to obtain the hit spectra of mesons such as 0, K and K® and proton 
fit function to obtain the mx spectra of baryons such as A, A, S and Q etc. 

All the mesons and baryons used in this analysis along with mode of measurement, 
rapidity and px ranges with their references are listed in Tabled for P+P and in Table HT1 for 



TABLE I. Particles with their mode of measurement, rapidity and px ranges for p+p collisions at 
different energies. 



Particle 


Mode 


Pt range (GeV) 


Rapidity range 


Reference 


p+p collision at >/s= 62.4 GeV 






77 
dE/dx 
dE/dx 


0.6-6.7 
0.3-2.85 
0.45-1.95 


\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.35 


PHENIX [17] 
PHENIX [16] 
PHENIX [16] 


p +p collision at y / s= 200 GeV 




P 

A 

n- 


77 
dE/dx 
dE/dx 
P ir- 
Air~ 
E- 7 
A if- 


0.6-18.9 
0.3-2.6 
0.468-6.5 
0.35-4.75 
0.67-3.37 

1.05-2.83 


\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.55 
\y\ < 0.75 
\y\ < 0.75 

\y\ < 0.75 


PHENIX [26] 
PHENIX [27] 
STAR ^8] 
STAR [191 
STAR [19] 

STAR ^9] 


p+p collision at y/s= 900 GeV 




K± 
P 

K° 
A 


dE/dx 
dE/dx 
dE/dx 

7T+7T — 
P 7T^ 
A7T~ 

s- 7 


0.1-2.5 
0.2-2.3 
0.3-2.3 
0.2-2.7 
0.7-3.25 
1.0-2.5 


\y\ < 0.8 
\y\ < 0.8 
|y| < 0.8 
\y\ < 0.8 
bl < 0.7 
\y\ < 0.8 


ALICE [24] 
ALICE [24] 
ALICE [24] 
ALICE [25] 
ALICE [25] 
ALICE [25] 



d+Au and Au+Au systems at different energies. The errors on the data are quadratic sums 
of statistical and uncorrelated systematic errors wherever available. Here we use all data from 
PHENIX (\y\ < 0.35), STAR (\y\ < 0.5, \y\ < 0.75) and ALICE (\y\ < 0.8) collaborations. 
The pion and proton spectra measured in p+p collisions for different energies are fitted using 
Eq. [1] and the parameters are given in Table II III and IIVI respectively. The parameters for 
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TABLE II. Particles with their measured decay channels, rapidity and px range for d+Au, Au+Au 
systems at y / i]v7V = 200 GeV. 



Particle 


Mode 


Pt range (GeV) 


Rapidity range 


Reference 


d+Au collision 




P 


dE/dx 


0.468-3.5 


\y\ 


< 0.75 


STAR [20] 


A++ 


N 7T 


0.3-1.5 


\y\ 


< 0.75 


STAR [21] 


Au +Au collision 




P 


dE/dx 


0.5-10.85 


\y\ 


< 0.75 


STAR [22J 


A 


P 71- 


0.65-4.75 


\y\ 


< 0.75 


STAP 1231 




E- 7 


0.85-4.75 


\y\ 


< 0.75 


STAR [23] 




A K~ 

H-7T 


1.31-4.0 


\y\ 


< 0.75 


STAR [23] 



TABLE III. The parameters of the Hagedorn distribution obtained by fitting pion spectra measured 
in p+p collisions at different center of mass energies. 



Parameters for 








pion 


= 62.4 GeV 


= 200 GeV [1JJ 


= 900 GeV 


A (GeV/c)- 2 


7.98 ± 0.42 


11.11 ± 0.37 


16.32 ± 0.81 


a (GeV/c)- 1 


0.137 ± 0.014 


0.32 ± 0.03 


0.82 ± 0.17 


b (GeV/c)" 1 


0.0 (fixed) 


0.024 ± 0.011 


0.0 (fixed) 


Po (GeV/c) 


1.24 ± 0.026 


0.72 ± 0.03 


0.43 ± 0.036 


n 


12.07 ± 0.11 


8.42 ± 0.12 


6.05 ± 0.28 



d+Au and Au+Au collisions are given in Table IVl and IVII respectively for pions and protons. 
For Au+Au system, data corresponding to four centrality classes namely 10-20 %, 20-40 %, 
40-60 % and 60-80 % has been analyzed. 
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TABLE IV. The parameters of the Hagedorn distribution obtained by fitting proton spectra mea- 
sured in p+p collisions at different center of mass energies. 



Parameters for 


V 


V 


V 


proton 


= 62.4 GeV 


= 200 GeV 


= 900 GeV 


A (GeV/c)" 2 


16.34 ± 5.847 


168.87 ± 45.77 


305.93 ± 38.69 


a (GeV/c)- 1 


0.51 ± 0.108 


0.15 ± 0.058 


0.283 ± 0.067 


b (GeV/c)- 1 


0.0 (fixed) 


0.0 (fixed) 


0.0 (fixed) 


Po (GeV/c) 


0.87 ± 0.104 


0.69 ± 0.029 


0.36 ± 0.007 


n 


11.49 ± 0.77 


9.86 ± 0.13 


7.17 ± 0.005 



TABLE V. The parameters of the Hagedorn distribution obtained by fitting pion spectra measured 
in d+Au and Au+Au (10-20%) collisions at y/s^ = 200 GeV. 



Parameters for 


d+Au system 


Au+Au system (10-20% cent.) 


pion 


Lii] 




A (GeV/c)" 2 


52.30 ± 1.63 


1564.36 ± 36.68 


a (GeV/c)- 1 


0.24 ± 0.01 


0.438 ± 0.007 


b (GeV/c)- 1 


0.11 ± 0.01 


0.213 ± 0.007 


Po (GeV/c) 


0.77 ± 0.02 


0.713 ± 0.004 


n 


8.46 ± 0.07 


8.383 ± 0.017 



TABLE VI. The parameters of the Hagedorn distribution obtained by fitting proton spectra mea- 



sured in d+Au (MB) and Au+Au collisions (for different centralities) at y/smq = 200 GeV. 



Parameters for proton 


d+Au collision 


Au+Au collision for different centralities 


10-20 % 


20-40 % 


40-60 % 


60-80 % 


A (GeV/c)" 2 


0.35 ± 0.10 


0.3 ± 0.2 


0.98 ± 0.13 


0.43 ± 0.13 


0.90 ± 0.34 


a (GeV/c)- 1 


1.63 ± 0.17 


1.07 ± 0.08 


0.93 ± 0.16 


1.12 ± 0.054 


1.16 ± 0.09 


b (GeV/c)- 2 


O.O(fixed) 


0.0 (fixed) 


O.O(fixed) 


O.O(fixed) 


O.O(fixed) 


Po (GeV/c) 


1.09 ± 0.045 


2.3 ± 0.12 


1.95 ± 0.18 


1.78 ± 0.051 


1.36 ± 0.05 


n 


7.06 ± 0.24 


10.22 ± 0.21 


10.34 ± 0.35 


9.48 ± 0.16 


9.05 ± 0.19 
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III. RESULTS AND DISCUSSIONS 



Figure (JT^i) shows the invariant yield of ^ 16 1, 7r° [l^ and (JTJd) shows the invariant yield 
of K ± 16J as a function of my measured in p+p collisions at y/s = 62.4 GeV. Here the data 
of neutral and charged pion are fitted with modified Hagedorn function shown in solid black 
curve. The corresponding fit parameters of pion are given in Table III IL In case of kaon, the 
solid line is obtained from pion fit function using mx scaling; the relative normalization has 
been used to fit the kaon spectra. There is difference in the measured data for K + and K~ 
and the tut scaled curve reproduces only K~ well. The relative normalization for kaons are 
given in Table IVHI . 

Figure |2] shows the invariant yield of different baryons as a function of rriT in p+p collisions 
at a/s = 200 GeV. Here the data of proton and anti-proton [3] are fitted with the modified 
Hagedorn function (shown in solid black curve). The fit parameters are given in Table [VTl 
which are used to get the ttit scaled spectra for measured baryons such as A, A, S _ , H+ and 



Q~ , Q + 19] as shown in figure. The relative normalizations of baryon (given in Table ITX]) to 
proton have been obtained by fitting the measured spectra. It is seen that the hit spectrum 
of all the baryons are well described by the m<r scaled curve. In case of A there are deviations 
of data from the curve in high region. 



The Figure ([3^) shows the invariant yield of proton and anti-proton 20( in d+Au system 
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m T - m [GeV/c 2 ] m T - m [GeV/c 2 ] 



? IG. 2. (Color online) The invariant yield of (a) proton (solid circle), anti-proton (open square) 
la], (b) A (solid circle), A (open square) [19], (c) H~ (solid circle), H+ (open square) [19| and (d) 



(fl- + J2+)/2 (solid circle) 



19(] as a function of my for p+p collision at yfs = 200 GeV. 



at v /J/vjv = 200 GeV as a function of m^. The solid black curve is fit function (Eq. [T]) the 
parameters of which are given in the Table I VI I The Figure (J3Jd) shows the invariant yield of 
A ++ [21] as a function of The solid line is tut scaled curve; the relative normalisation 
given in Table HX] is obtained by fitting the measured data reproducing the shape of the 
measured spectrum very well. 



Figure H] shows the invariant yield of proton and anti-proton 20] as a function of % for 



Au+Au collsion at y / J/v r /v = 200 GeV for different centralities. The solid black curves are 
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d+Au at ^fs^ = 200 GeV 
• P 

□ p 

— Fitted curve 




0.5 



1 1.5 2 2.5 
m T - m [GeV/c 2 ] 



10- 1 



> 
o 



Q. 



d+Au at \/s^ = 200 GeV 
• A"(1232) 
— m T scaled curve 

Fitted norm. = 0.497 ± 0.007 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
m T - m [GeV/c 2 ] 



IG. 3. (Color online) The invariant yield of (a) proton (solid circle), anti-proton (open square) 



20J and (b) A ++ baryons 



2l| as a function of niT for d+Au collision at ^/snn = 200 GeV. 




FIG. 4. (Color online) The invariant yield of proton (closed symbols) and anti-proton (open 
symbols) [22] as a function of mj- in Au+Au collisions at ^snn = 200 GeV for different centralities. 
The solid lines are Hagedorn fit function. 
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Au+Au at \/s NN = 200 GeV; 
A (closed symbols) 

A (open symbols) - 

10-20% centrality : 

20-40% centrality ~| 
40-60% centrality 

60-80% centrality — 

m T scaled curves = 




1 1.5 2 2.5 3 

ny m [GeV/c 2 ] 



FIG. 5. (Color online) The invariant yield of A (closed symbols) and A (open symbols) 23|] as 
a function of mx in Au+Au collisions at y / J/vjv = 200 GeV for different centralities. The solid 
lines are obtained using scaling; the relative normalization has been used to fit the measured 
spectra. 



the fitted modified Hagedorn functions with parameters given in Table I VI I Figure shows 
the invariant yield of A and A [23] for different centralities as a function of % for Au+Au 
collision at a/saT/v = 200 GeV. The solid black curves are the m T scaled curves obtained 
from the proton spectra of corresponding centralities with a fitted normalization given in 
Table IIXI Figure [6] shows the invariant yield of H and S [23( for different centralities as a 
function of vtlt for Au+Au collision at ^/snn = 200 GeV. The solid black curves are the nix 
scaled curves obtained from the proton spectra of corresponding centralities with a fitted 
normalization given in Table HXl Here we see that, in both the cases for A, A and E+, the 
experimental data are in good agreement with the % scaled curve for the most peripheral 
centrality bin. But for other centralities, the nature of the m-r scaled curves are different 
from the shape of data. Same obeservation is made from Figure [7] which shows the invariant 
yield of f2 |23| as a function of tut for Au+Au collision at ^/s^n = 200 GeV. 

Figure (jHfci) shows the invariant yield of 7r ± |24| as a function of m T for p+p collision at 
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Au+Au at = 200 GeV 
S" (closed symbols) 
3" (open symbols) e 
10-20% centrality ; 
20-40% centrality -= 
40-60% centrality = 
60-80% centrality 
scaled curves 




1 1.5 2 2.5 3 

ny m [GeV/c 2 ] 



FIG. 6. (Color online) The invariant yield of H~ (closed symbols) and H+ (open symbols) 23J] as 
a function of mj- in Au+Au collisions at ^/sJjn = 200 GeV for different centralities. The solid 
lines are obtained using scaling; the relative normalization has been used to fit the measured 
spectra. 



y/s = 900 GeV. The solid black curve is fit function (Eq. [T]) whose parameters are given in 
the Table UTTl Figures [8] (b), (c) and (d) show the invariant yields of 24(, Kg 25|, and 



25j, respectively as a function of mj. The solid black curves are the ttit scaled curves 



obtained from the pion spectrum with a fitted normalization given in Table IVHI From this 
we clearly see that the wit scaled curves do not reproduce mesons for p+p at 900 GeV that 
well. 



The Figure M shows the invariant yields of proton and anti-proton 24| as a function of m,T 
for p+p collision at yfs = 900 GeV. The solid black curve is the fitted modified Hagedorn 
function with parameters given in Table IIVI The Figures |9] (b) and (c) show the invariant 
yields of A [25j and S 25|. The solid black curves are the % scaled curves obtained from 



the proton spectra with fitted normalizations given in the Tabe IVIIII Here the scaled 
curve is in good agreement with A and with small data set of H at \/s = 900 GeV. 

Figure [10] shows the invariant yields of proton (P), anti-proton (P) and pion (jr) as a 
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Au+Au at xfs^ = 200 GeV 
£2 — > A K" 
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FIG. 7. (Color online) The invariant yield of O (closed symbols) [23(] as a function of mx in 
Au+Au collisions at ^s^n = 200 GeV for different centralities. The solid lines are obtained using 
wit scaling; The Relative normalization has been used to fit the measured spectra. 



function of m T for p+p collisions at y/s = 62.4, 200 and 900 GeV and for Au+Au collisions 
(10-20% centrality) at ^snn = 200 GeV along with their fit functions obtained from Eq. [TJ 
The shapes of pions and protons m-r spectra are different for all systems at all energies. In 
comparision to the proton and pion production at low energies (e.g. 62.4 and 200 GeV), the 
proton production is less as compared to the pion production at 900 GeV. 



Figure [TT] shows the ratio of proton to pion (p/7r) (obtained from the fit functions in 
Figure [TO]) vs m T for p+p collisions at y/s = 62.4, 200 and 900 GeV and for Au+Au 
collisions (10-20% centrality) at ^/snn = 200 GeV. It is seen that for p+p collisions, the 
shape of p/n curves are same at 62.4, 200 and also at 900 GeV, but the ratio of p/V is very 
small at ^/snn = 900 GeV. The enhancement of p/-7r ratio at intermediate m<r in Au+Au 



is understood in terms of recombination model 
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p+p at\fs= 900 GeV 

— m T scaled curve 

Fitted norm. = 0.445 + 0.006 
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— m T scaled curve 

Fitted norm. = 0.230 ± 0.019 
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FIG. 8. (Co 
solid circle) 



or online) The invariant yield of (a) ir^ (closed triangle) 24], (b) K.^ (open square, 



24j, (c) Kg (solid circle) [25]] and (d) <p (solid circle) [25j as a function of itit for p+p 



collision at yfs = 900 GeV. The solid lines are obtained using tut scaling. 



TABLE VII. The relative normalization (S) of mesons obtained by fitting the my scaled spectra 
with the measured spectra for p+p collisions at -y/s = 62.4, 200 and 900 GeV. 



s 


y/i = 62.4 GeV 


y/s = 200 GeV 


^f~s = 900 GeV 






From [11] 




K/tt 


0.38 ± 0.02 


0.422 ± 0.014 


0.424 ± 0.02 








0.44 ± 0.02 


0/7T 




0.233 ± 0.013 


0.23 ± 0.02 
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m T -m [GeV/c 2 ] m T -m [GeV/c 2 ] 




iTG. 9. (Color online) The invariant yield of (a) proton (solid circle), anti-proton (open square) 
24|, (b) A (solid circle), A (open square) [25j and (c) S (solid circle) [25j as a function of my for 
p+p collision at y/s = 900 GeV. The solid lines are obtained using mx scaling. 



TABLE VIII. The relative normalization (S) of baryons obtained by fitting the m-r scaled spectra 
with the measured spectra for p+p collisions at different center of mass energies. 



s 








= 200 GeV 


= 900 GeV 


A/p 


0.785 ± 0.021 


0.93 ± 0.2 


H/P 


0.14 ± 0.03 


0.2 ± 0.02 




0.060 ± 0.022 
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FIG. 10. (Color online) The invariant yields of proton, anti-proton, tt^ and ir° as a function of 
m T for (a) p+p at ^fs = 62.4 GeV Q, Q, (b) p+p at ^ = 200 GeV Q, fl 27 L (c) p+p at yfs 



900 GeV 



2J, (d) Au+Au (for 10-20% centrality) at y/sjfjf = 200 GeV [2i 



28|, 



29J. The solid 



lines show the modified Hagedorn fit function drawn over proton and pion. 



IV. CONCLUSION 



In summary, we have made a systematic study of transverse mass spectra of mesons and 
baryons at RHIC and LHC energies. For p+p and d+Au collision at ^snn = 200 GeV the 
baryon rrix spectra scale with proton tut spectrum just like all meson spectra scale with 
pion spectra. In p+p 200 GeV, ttit scaling reproduce data well for all the baryons except 
A as there is a disagreement at high points. In case of Au+Au collisions, the strange 
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For p+p collision 




m T [GeV/c 2 ] 

FIG. 11. (Color online) The ratio of proton fit function to the pion fit function (p/vr) as a function 
of rriT for p+p collisions at different energies and Au+Au collisions (10-20% centrality) at 200 
GeV. 



TABLE IX. The relative normalization (S) of baryons obtained by fitting the tut scaled spectra 



with the measured spectra for different collision systems at ^/snn = 200 GeV. 



s 


p+p collision 


d+Au collision 


Au+Au collision for different centralities 


10-20% 


20-40 % 


40-60 % 


60-80 % 


A/P 


0.785 ± 0.021 




0.70 ± 0.09 


0.70 ± 0.09 


0.7 ± 0.2 


0.75 ± 0.28 


S/P 


0.140 ± 0.030 




0.18 ± 0.02 


0.18 ± 0.02 


0.18 ± 0.05 


0.17 ± 0.06 


fi/p 


0.060 ± 0.022 




0.16 ± 0.008(0-5%) 


0.13 ± 0.01 


0.12 ± 0.03 




A/p 




0.497 ± 0.016 











baryons behave differently from protons for all centralities but there is better agreement in 
the two for peripheral collisions. This study has also been performed on the meson and 
baryon spectra in p+p collisions at 62.4 GeV which lies in between the highest RHIC energy 
and SPS energy where the scaling was first observed and we find that the ttit scaling 
of mesons work at 62.4 GeV (Baryon data at 62.4 GeV is not available). We test the 
scaling at LHC with the first meson and baryon spectra measured in p+p collisions at 900 
GeV and find that the uit spectra of kaons and 4> do not scale with pions that well. At 
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900 GeV, A is produced well by my scaling where as the S spectrum is reproduced well 
within the limited availability of data points. The shapes of pions and protons nix spectra 
are different at all energies. Proton to pion ratio is studied at different energies (using fit 
functions) and for different systems and is found that for p+p collisions, the behaviour of 
p/n curves are similar at 62.4, 200 and also at 900 GeV. It will be interesting to have such 
study at higher energies such as 2.76 TeV and 7 TeV, once the data tables are available. 
More comprehensive measurements at lower energies such as 62.4 GeV will benifit to draw 
a systematic picture of hadron production as a function of energy. 
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